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Abstract 

The time evolution equations of a simplified isolated ideal gas, the "tetra- 
hedral" gas, are derived. The dynamical behavior of LMC (Lopez, Mancini, 
Calbet) complexity is studied in this system. In general, it is shown that the 
complexity remains within the bounds of minimum and maximum complex- 
ity. We find that there are certain restrictions when the isolated "tetrahedral" 
gas evolves towards equilibrium. In addition to the well-known increase in en- 
tropy, the quantity called disequilibrium decreases monotonically with time. 
Furthermore, the trajectories of the system in phase space approach the max- 
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imum complexity path as it evolves toward equilibrium. 



Typeset using REVT^ 



I. INTRODUCTION 



Several definitions of complexity, in the general sense of the term, have been presented 
in the literature. These can be classified according to their calculation procedure into two 
broad and loosely defined groups. 

One of these groups is based on computational science and consists of all definitions 
based on algorithms or automata to derive the complexity. Examples are the logical depth 
||T|], the e- machine complexity and algorithmic complexity [0. These definitions have 
been shown to be very useful in describing symbolic dynamics of chaotic maps, but they 
have the disadvantage of being very difficult to calculate. 

Another broad group consists of those complexities based on the measure of entropy or 
entropy rate. Among these, we may cite the metric or K-S entropy rate the thermody- 
namic depth 1^ , the effective measure complexity [0 , and the simple measure for complexity 
0. These definitions have also been very useful in describing symbolic dynamical maps, the 
latter having been applied to a non-equilibrium Fermi gas 0. They suffer the disadvantage 
of either being very difficult to calculate or having a simple relation to the regular entropy. 

New definition types of complexity have recently been introduced. These are based 
on quantities that can be calculated directly from the distribution function describing the 
system. One of these is based on "meta-statistics" |]TU| and the other on the notion of "dise- 
quilibrium" |]TI| . This latter definition will be referred to hereafter as the LMC complexity. 
These definitions, together with the simple measure for complexity described above, have 
the great advantage of allowing easy calculations within the context of kinetic theory and 
of permitting their evaluation in a natural way in terms of statistical mechanics. 

The disequilibrium-based complexity is easy to calculate and shows some interesting 
properties [|lT|, but suffers from the main drawback of not being very well behaved as the 



system size increases, or equivalently, as the distribution function becomes continuous [12 



Feldman and Crutchfield tried to solve this problem by defining another equivalent term for 
disequilibrium, but ended up with a complexity that was a trivial function of the entropy. 
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Whether these definitions of complexity are useful in non-equilibrium thermodynamics 
will depend on how they behave as a function of time. There is a general belief that, although 
the second law of thermodynamics requires average entropy (or disorder) to increase, this 



does not in any way forbid local order from arising |T3|]. The clearest example is seen 
with life, which can continue to exist and grow in an isolated system for as long as internal 
resources last. In other words, in an isolated system the entropy must increase, but it should 
be possible, under certain circumstances, for the complexity to increase. 

In this paper we will examine how LMC complexity evolves with time in an isolated 
system and we will show that it indeed has some interesting properties. The disequilibrium- 



based complexity defined in Ref. [|TT| actually tends to be maximal as the entropy increases 
in a Boltzmann integro-differential equation for a simplified gas. 

In Sec. H LMC complexity definition is reviewed. We proceed to calculate the distribu- 
tions which maximize and minimize the complexity and its asymptotic behavior, and also 
introduce the basic concepts underlying the time evolution of LMC complexity in Sec. fT|. 
Later, in Sec. fV|, by means of numerical computations following a restricted version of the 
Boltzmann equation, we apply this to a special system, which we shall term "tetrahedral 
gas". Finally, in Sec. [V|, the results and possible future lines of investigation are discussed, 
together with their possible applications. In the Appendix is given an analytical and a 
numerical demonstrations of the results of the numerical calculations for the tetrahedral 
gas. 



II. DEFINITION OF COMPLEXITY 

The definition of LMC complexity, C, is 

C = D-H, (1) 
where D is the disequilibrium term and H is the entropy. 

We assume that the system can be in one of possible accessible states, i. The prob- 
ability of the system being in state i will be given by the discrete distribution function, /j. 
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The system is defined such that, if isolated, it will reach equilibrium, with all the states 
having equal probability, f^. 

These definitions imply that all values of fi are positive, 

fi > 0, (2) 
and that a normalization, /, must hold such that 

iV 

^ = E/^ = 1' (3) 

i=l 

and the equilibrium distribution function is, 

= ^- w 

The definition of disequilibrium, D, is given as a distance to the probability in equilib- 
rium, /e: 

N 



i=l 

The normalized entropy, H, is defined as 

1 ^ 

H = y/iln/i. (6) 

Note that since < < 1 and < D < {N — 1)/N, the complexity, C, is normalized 
(0<C< 1). 

III. COMPLEXITY VERSUS TIME 

A. Complexity versus entropy 

We are interested in an isolated system with an initial arbitrary discrete distribution, 
and which evolves toward equilibrium, in which it reaches an equiprobability distribution. 

To study the time evolution of the complexity, a diagram of C versus time, t, can be 
used. But, as we know, the second law of thermodynamics states that the entropy grows 
monotonically with time; that is. 



This implies that an equivalent way to study the time evolution of the complexity can 
be obtained by plotting C versus H. In this way, the entropy substitutes the time axis, 
since the former increases monotonically with the latter. The conversion from C vs. H to 
C vs. t diagrams is achieved by stretching or shrinking the entropy axis according to its 
time evolution. This method is a key point in all this discussion. Note that, in any case. 



the relationship of H versus t will, in general, not be a simple one |14 



B. Maximum and minimum complexity 

When an isolated system evolves with time, the complexity cannot have any possible 
value in a C versus H map, but it must stay within certain bounds, Cmax and Cmin- These 
are the maximum and minimum values of C for a given H. Since C = D ■ H, finding the 
extrema of C for constant H is equivalent to finding the extrema of D. 

There are two restrictions on D: the normalization, J, and the fixed value of the entropy, 
H. To find these extrema undetermined Lagrange multipliers are used. Differentiating Eqs. 
(|), (H) and d^) we obtain, 

dD 

^ = 2(/,-/e), (8) 



dH _ 1 
'dfj ~ ~lniV 



— = -TrT7(in/, + i). (10) 



Defining Ai and A2 as the Lagrange multipliers, we get: 

2{fj - /e) + Ai + A2(ln/, + 1)/ IniV = 0. (11) 

Two new parameters, a and f3, which are a linear combinations of the Lagrange multi- 
pliers are defined: 

fj + a\nf,+(3 = 0, (12) 



where the solutions of this equation, fj, are the values that minimize or maximize the 
disequilibrium. 

In the maximum complexity case there are two solutions, fj, to Eq. (|1^) which are shown 
in Table |. One of these solutions, /max, is given by 

1 



H 



IniV 



/max ll^ /max ~l~ /max) 1-^1 



N-1 



(13) 



and the other solution by (1 — fma.x)/iN — 1). 

The maximum disequilibrium, -Dmax, for a fixed H is 

/^max=(/max-/e)' + (iV-l)(^^^-/e) , (M) 

and thus, the maximum complexity, which depends only on H, is 

Cmax(-f^) = -Dmax " H. (15) 

Note that the behavior of the maximum value of complexity versus In N has been studied 
in Ref. 

Equivalently, the values, fj, that give a minimum complexity are shown in Table One 
of the solutions, /min, is given by 

1 



H 



InN 



/min 111 /min ~l~ (1 /min) 1^1 



N -n-1 



(16) 



where n is the number of states with fj = and takes a value in the range n = 0, 1, ... , N—2. 
The resulting minimum disequilibrium, -Dmin, for a given H is, 

I^mm = (/min - fe)' + (N - U - 1) (^^^^^ " f^j + ^/e • (17) 

Note that in this case fj = is an additional hidden solution that stems from the positive 
restriction in the /j values, Eq. (0). To obtain these solutions explicitly we can define Xi 
such that 

f^ = (18) 
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These Xi values do not have the restriction imposed by Eq. and can take a positive 
or negative value. If we repeat the Lagrange multiplier method with these new variables a 
new solution arises: xj = 0, or equivalently, fj = 0. 

The resulting minimum complexity, which again only depends on H, is 

CminiH) = Anin ■ H. (19) 

As an example, the maximum and minimum of complexity, C^^x and Cmin, are plotted 
as a function of the entropy, H, in Fig. ^ for = 4. In Fig. ^ the maximum and minimum 
disequilibrium, Z^max and -Dmin, versus H are also shown. 



C. Minimum "envelope" 

As shown in Fig. |^ the minimum disequilibrium function is piecewise defined, having 
several points where its derivative is discontinuous. Each of these function pieces corresponds 
to a different value of n (Table |I|) . In some circumstances it might be helpful to work with the 
"envelope" of the minimum disequilibrium function. The function, -Dminenv, that traverses 
all the discontinuous derivative points in the i^min versus H plot is 

D ■ - f-^^""^ - — (20) 

-'-^minenv — V 

and is also shown in Figure 0. 

D. Asymptotic behavior of the complexity as N — ^ oc 

When A^ tends toward infinity the probability, /max, of the dominant state has a linear 
dependence with the entropy, 

lim /„,ax =l-H, (21) 

N^oo 

and thus the maximum disequilibrium scales as 

lim D^,, = (1 - Hf. (22) 



The maximum complexity tends to 



lim C^a. = H-{1- Hf. (23) 



The hmit of the minimum disequihbrium and complexity vanishes, 



lim Dmincnv = 0, (24) 
Af^oo 



and thus 



lim Cmin = 0. (25) 

Af— >oo 

In general, in the limit N ^ oo, the complexity is not a trivial function of the entropy, in 
the sense that for a given H there exists a range of complexities between and Cmax (Eqs. 
m and H). 

In particular, in this asymptotic limit, the maximum of Cmax is found when H = 1/3, 
or equivalently /max = 2/3, which gives a maximum of the maximum complexity of Cmax = 
4/27. This was numerically calculated by Anteneodo and Plastino in Ref. |TH. 



IV. AN EXAMPLE: THE TETRAHEDRAL GAS 
A. The tetrahedral gas 

We present a simplified example of an ideal gas: the tetrahedral gas. This system is 
generated by a simplification of the Boltzmann integro-differential equation of an ideal gas. 
We are interested in studying the disequilibrium time evolution. 

The Boltzmann integro-differential equation of an ideal gas with no external forces and 
no spatial gradients is 
a/(v;t) 



dt 



J J dQ,,^a{v, - V ^ v: - v') I V, - v| [/(v;; t)/(v'; t) - /(v,; t)/(v; t)] , 

(26) 
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where a represents the cross section of a colhsion between two particles with initial velocities 
V and V* and after the collision with velocities v' and v'^; and fic.m. are all the possible 
dispersion angles of the collision as seen from its center of mass. 

In the tetrahedral gas, the particles can travel only in four directions in three-dimensional 
space and all have the same absolute velocity. These directions are the ones given by joining 
the center of a tetrahedron with its corners. The directions can be easily viewed by recalling 
the directions given by a methane molecule, or equivalently, by a caltrop, which is a device 
with four metal points so arranged that when any three are on the ground the fourth projects 
upward as a hazard to the hooves of horses or to pneumatic tires (see Fig. ^). 

By definition, the angle that one direction forms with any other is the same. It can be 
shown that the angles between different directions, a, satisfy the relationship cos a = —1/3, 
which gives a = 109.47°. The plane formed by any two directions is perpendicular to the 
plane formed by the remaining two directions. 

We assume that the cross- sect ion, a, is different from zero only when the angle between 
the velocities of the colliding particles is 109.47°. It is also assumed that this collision 
makes the two particles leave in the remaining two directions, thus again forming an angle 
of 109.47°. A consequence of these restrictions is that the modulus of the velocity is always 
the same no matter how many collisions a particle has undergone and they always stay 
within the directions of the vertices of the tetrahedron. Furthermore, this type of gas does 
not break any law of physics and is perfectly valid, although hypothetical. 

We label the four directions originating from the center of the caltrop with numbers, 
1, 2, 3, 4 (see Fig. |]). The velocity components with the same direction but opposite sense, 
or equivalently, directed toward the center of the caltrop, are labeled with negative numbers 
-1,-2,-3,-4. 

In order to formulate the Boltzmann equation for the tetrahedral gas, and because all 
directions are equivalent, we need only study the different collisions that a particle with one 
fixed direction can undergo. In particular if we take a particle with direction —1 the result 
of the collision with another particle with direction —2 are the same two particles traveling 
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in directions 3 and 4; that is, 

(-1,-2) ^ (3,4). 
With this in mind the last bracket of Eq. ( PB| ) is, 

73/4 — /-1/-25 

where /j denotes the probabihty of finding a particle in direction i. Note that the dependence 
on velocity, v, of the continuous velocity distribution function, /(v; t), of Eq. (^6]) is in our 
case contained in the discrete sub index, i, of the distribution function /j. 
We can proceed in the same manner with the other remaining collisions, 

(-1,-3)^(2,4), 
(-1,-4) ^(2,3). 

When a particle with direction —1 collides with a particle with direction 2, they do not 
form an angle of 109.47°; i.e., they do not collide, they just pass by each other. This is 
a consequence of the previous assumption for the tetrahedral gas, which establishes a null 
cross section for angles different from 109.47°. The same can be said for collisions (—1,3), 
(—1, 4), and (—1, 1). AH these results are summarized in Table |lll. 



Taking all this into account, Eq. (^6]) for direction —1 is reduced to a discrete sum, 

"^^"i = (/3/4 - /-1/-2) + (/2/4 - /-1/-3) + (/2/3 - /~l/-4), (27) 



dt 

where all other factors have been set to unity for simplicity. 
The seven remaining equations are: 

df-2 



dt 

df-3 

dt 
df-4 

dt 



ifsh - /-1/-2) + (/1/4 - /-2/-3) + (/1/3 - /-2/-4), 
(/2/4 - /-3/-1) + (/4/1 - /-3/-2) + (/1/2 - /-3/-4), 
(/2/3 - /-4/-1) + (/3/1 - /-4/-2) + (/1/2 - /-4/-3), 



^ = (/_3/-4 - /1/2) + (/-2/-4 - /1/3) + (/-2/-3 - /1/4) 
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^ = (/-2/-4 - /3/1) + - /3/2) + if-lU - hh), 

^ = (/-2/-3 - /4/1) + (/-1/-3 - /4/2) + (/-1/-2 - /3/4). (28) 

If we now make fi = f-i{i = 1,2,3,4) initially, this property is conserved in time. The 
final four equations defining the evolution of the system are: 

^ = (/3/4 - /1/2) + (/2/4 - /1/3) + (/2/3 - /1/4), 
^ = (/3/4 - /1/2) + (/1/4 - /2/3) + (/1/3 - /2/4), 
f H/2/4-/3A).(A/.-/3/2) + (A/2-/3/4), 

^ = (/2/3 - /4/1) + (/1/3 - /4/2) + (/1/2 - /3/4). (29) 

Note that the ideal gas has been reduced to the tetrahedral gas, which is a four- 
dimensional dynamical system. The velocity distribution function, corresponds to the 
probability distribution function of Sec. |I| with = 4 accessible states. 



B. Evolution of the tetrahedral gas with time 



The tetrahedral gas (Eqs. p9[) reaches equilibrium when fi = 1/N for z = 1, 2, 3, 4 and 
= 4. This stationary state, dfi/dt = 0, represents the equiprobability towards which the 

system evolves in time. This is consistent with the previous definition of disequilibrium, Eq. 

(^), in which we assumed that equilibrium was reached at equiprobability, where D = 0. 
As the isolated system evolves it gets closer and closer to equilibrium. In this sense, 

one may intuitively think that the disequilibrium will decrease with time. In fact, it can 

be shown that, as the system approaches equilibrium, D tends to zero monotonically with 

time: 

dD , , 

The analytical demonstration of this inequality for the tetrahedral gas is shown in Ap- 
pendix 1^. 
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There are even more restrictions on the evolution of this system. It would be expected 
that the system approaches equilibrium, D = 0, hj following the most direct path. To 
verify this, numerical simulations for several initial conditions have been undertaken. In all 
of these we observe the additional restriction that D approaches -Dmax on its way to D = 0. 
In fact it appears as an exponential decay of D towards -Dmax in a D versus H plot. As an 
example, two of these are shown in Fig. ^, where Fig. ^(a) shows a really strong tendency 
towards -Dmax- Contrary to intuition, among all the possible paths that the system can 
follow toward equilibrium, it chooses those closest to -Dmax in particular. 

We can also observe this effect in a complexity, C, versus H plot. This is shown for the 
same two initial conditions as the previous figure in Figure |^. 

This additional restriction to the evolution of the system is better viewed by plotting 
the difference Cmax — C versus H. In all the cases analyzed(see two in Fig. ||) the following 
condition is observed: 

<^(Cmax ~ C) 



dt 



< 0. (31) 



This has been verified numerically and is illustrated in Figure |^, where this time deriva- 
tive, which always remains negative, is shown as a function of H for a grid of uniformly 
spaced distribution functions, (/i, /2, /a, /a), satisfying the normalization condition, Eq. (H). 
Two system trajectories are also shown for illustrative purposes. The numerical method 
used to plot this function is explained in Appendix 0. 

C. Maximum complexity path as an attractive trajectory 

As shown in Table |I[ a collection of maximum complexity distributions for = 4 can 
take the form 



fl /max 

/, = i^i^,^ = 2,3,4 (32) 
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where /max runs from (equiprobability distribution) to 1 ("crystal" distribution). The 
complexity of this collection of distributions covers all possible values of Cmax- 

There is actually a time evolution of the tetrahedral gas, or trajectory of the system, 
formed by this collection of distributions. Inserting Eqs. ( |32D in the evolution Eqs. (p9D, it 
is found that all equations are compatible with each other and the dynamical equations are 
reduced to the relation, 

^ = ^(4/lx-5/max + l). (33) 

This trajectory is denoted as the maximum complexity path. 

Note that the equiprobability or equilibrium, /max = 1/4, is a stable fixed point and 
the maximum disequilibrium "crystal" distribution, /max = 1, is an unstable fixed point. 
Thus the maximum complexity path is a heteroclinic connection between the "crystal" and 
equiprobability distributions. 

The maximum complexity path is locally attractive. Let us assume, for instance, the 
following perturbed trajectory 

/l /max; 

1 /max 



/^ 



2 



3 



f3 = ^-^ + S, 

h = ^—j^ - (34) 



whose evolution according to Eqs. ( pof ) gives the exponential decay of the perturbation, 5: 

d6 /4/max + 2^ 



dt V 3 
showing the attractive nature of these trajectories. 



6, (35) 



V. CONCLUSION 



The time evolution of the LMC complexity, C, has been studied for a simplified model 
of an isolated ideal gas: the tetrahedral gas. In general, the dynamical behavior of this 
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quantity is bounded between two extremum curves, Cmax and Cmin, when observed in a C 
versus H phase space. These complexity bounds have been derived. 

For the isolated tetrahedral gas two constraints on its dynamics are found. The first, 
which is analytically demonstrated, is that the disequilibrium, D, decreases monotonically 
with time until it reaches the value D = for the equilibrium state. The second is that 
the maximum complexity paths, Cmax, are attractive in phase space. In other words, the 
complexity of the system tends to equilibrium always approaching those paths. We verify 
this numerically; that is, the time derivative of the difference between Cmax and C is negative. 

Fig. ^ summarizes the dynamical behavior of the tetrahedral gas. The different trajec- 
tories starting with arbitrary initial conditions, which represent systems out of equilibrium, 
evolve towards equilibrium approaching the maximum complexity path. 

Whether these properties are useful in real physical systems will depend on further work 
on this subject. More can be said about the macroscopical nature of the disequilibrium 
when this work is extended to more general systems, such as to the ideal gas following 
the complete Boltzmann integro-differential equation. Another feature, which could result 
useful, would be to approximate the evolution of a real physical system trajectory to its 
maximum complexity path. Note that in general, for a real system, the calculation of the 
maximum complexity path will not be an easy task. 

APPENDIX A: PROOF OF THE MONOTONICALLY DECREASE IN 
DISEQUILIBRIUM WITH TIME 



We now present the analytical proof of the inequality shown in Eq. (^) from Sec. ^VB 

The time derivative of the disequilibrium, for the isolated tetrahedral gas, is explicitly 
j;iven by. 
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Using the normahzation Eq. we find that 

= 0, 

. dt 



and we are left with 

dt dt' 

If we substitute Boltzmann Eq. (pQ]) in the previous one, we obtain: 
F = 2~dt ^ ^^^'^^^ " •^i'^^ ~'~ •^i'^2/4 — /1/3 + /1/2/3 — /1/4 

+ /2/3/4 — /l/l + /1/2/4 — /s/l + /1/2/3 — 71/4 
+ /1/3/4 — 73/2 + 72/3/4 — /1/3 + /1/2/3 — /3/4 
+ /1/2/4 — /s/l + /1/3/4 — /2/4 + /2/3/4 — /l/l) 

where the new variable F has been defined. 

We now split this function F into two different terms: 

- flh - /l/l - /l/l - /I/2 

— /4/1 ~" /4/2 
+ /1/3/4 + 2/1/2/3 
+ 2/1/2/4 + /2/3/4 



and 



— flu - /I/3 

— /l /3 — /l /4 

— /l/s — /I/4 

+ /1/2/3 + 2/1/3/4 
+ 2/2/3/4 + /1/2/4. 



Then 
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F = G + K. (A8) 

Note that G and K are symmetrical functions, in the sense that we can transform one 
into the other by changing (/i, /a, /s, fi) by (/g, /4, /i, /a), respectively. 
To prove that G < 0, the following variable change is performed: 

/i = yi, 

/s = Z/3 + Z/4, 

fA = y3- va. (A9) 

The positivity of the distributions functions, Eq. (||), implies that, yi > and 1/2 > 0. 
The function G in the new variables reads, 

G = -yly2 - yhi - ylvi - ^ylvi - yly2 - ^yly2 + ^ymy^- (Aio) 

Regrouping terms, G can be expressed as 

G = -y2iyi - yzf - y\{y2 - yzf - ^ylyi - ^yjyi- (All) 

Since yi and 1/2 are both positive and the squared quantities are also positive, we conclude 
that 

G < 0. (A12) 

The same inequality can be demonstrated for K due to its symmetry with G, thus proving 
the assumption from Eq. ( [Al|) . 



APPENDIX B: METHOD OF CALCULATING THE TIME DERIVATIVE OF 
THE MAXIMUM COMPLEXITY MINUS THE COMPLEXITY 

To calculate the quantity 



17 



dt ^^^^ 

from some given values of the distribution functions, (/i, /2, /s, /4), we derive the expression 
of maximum complexity, Eq. minus the definition of complexity, Eq. ([^), obtaining, 

d{C^..-C) dD^, dH dD dH 

dt = -IT^ + ^--^ -IF^- 

Let us now examine each of these terms beginning from the end. The time derivative of 
the entropy is calculated by differentiating its definition, Eq. (P): 



dH 1 ^ 



dt IniV 
But recalling Eq. ( [A 3D we can simplify this to 

'df, 



dfi. r- . dfi 



(B3) 



dH 1 ^ 



(B4) 



dt IniV ^ ^. 

This last term can be easily calculated using the evolution Eqs. (pQ]). 

The terms D and H can be readily calculated using their definitions from Eqs. and 

(!)■ 

The time derivative of the disequilibrium has already been expressed in Appendix ^ with 
Eq. (^). 

The maximum disequilibrium is calculated with the previous found value of H by in- 
verting Eq. (|13|) and obtaining /max- After this, the expression defining the maximum 
disequilibrium, Eq. (0), can be used. 

Finally the time derivative of the maximum disequilibrium is calculated as follows: 

dDraax dDjj^g^ df^g^ dH (B5) 

dt C?/max dH dt 

The first two derivatives on the right-hand side can be calculated analytically with Eqs. 
(|1^) and (]T3|). The numerical value can be found using the previously calculated value of 

/max- 
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FIGURES 

FIG. 1. Maximum, minimum, and minimum envelope complexity, Cmax; C'min) and Cminenv 
respectively, as a function of the entropy, H, for a system with A*" = 4 accessible states. 

FIG. 2. Maximum, minimum, and minimum envelope disequilibrium, Dmax; -Cmim and -Dminenv 
respectively, as a function of the entropy, H, for a system with = 4 accessible states. 

FIG. 3. The four possible directions of the velocities of the tetrahedral gas in space. Positive 
senses are defined as emerging from the center point and with integer numbers, 1,2,3,4. 

FIG. 4. Time evolution of the system in (H, D) phase space for two different initial conditions 
at time t = 0: (a) (/i, /a, /s, /4) = (0.8,0.2,0,0) and (b) ifi,f2,h,h) = (0.5,0.5,0,0). The 
maximum and minimum disequilibrium are shown by dashed lines. 

FIG. 5. Time evolution of the system in (H, C) phase space for two different initial conditions 
at time t = 0: (a) (/i, /a, /s, /4) = (0.8,0.2,0,0) and (b) (/i, /a, /s, /4) = (0.5,0.5,0,0). The 
maximum and minimum complexity are shown by dashed lines. 

FIG. 6. Time evolution of the system in (/?, Cmax ~ C) phase space for two different initial 
conditions at time t = 0: (a) (fi, 12,13,14) = (0.8,0.2,0,0) and (b) (fi, f2, fs, Ia) = (0.5,0.5,0,0). 
The values Cmax — C'min are shown by dashed lines. 

FIG. 7. Numerical verification of (i(Cinax — C)/dt < 0. This time derivative is shown as a 
function of H. A grid of uniformly spaced, A/j = 0.5, distribution functions, {fi, f2, fs, Ia), 
satisfying the normalization condition Eq. (Q), have been used. Two system trajectories for initial 
conditions, t = 0, {fi, f2, fs, fi) = (0.8,0.2,0,0) and {fi, f2, fs, fi) = (0.5,0.5,0,0) are also shown 
for illustrative purposes. It can be seen how the above-mentioned time derivative always remains 
negative. 
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FIG. 8. Summary of this paper. The time evolution of the system for three different initial 
conditions, t = 0, (/i, /2, /s, /a) = (0.8, 0.2, 0, 0), (/i, /2, fs, fi) = (0.5, 0.5, 0, 0), and the maximum 
complexity path are shown. The minimum complexity is shown by dashed lines. It can be seen 
how the system tends to approach the maximum complexity path as it evolves in time toward 
equilibrium. 
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TABLES 

TABLE I. Probability values, fj. thai, givo a maximum of disof|uilibrium. -Duiax- for a givou II 



Number of states with 




fj 


Range of fj 


with fj 








1 

A^- 1 




/max 
1 /max 

A^- 1 


A^ 

1 


TABLE IL Probabihty values, fj, 


that give a minimum of disequilibrium, -Dmin, for a given H 


Number of states 




fj 


Range of fj 


with fj 








n 










1 

N-n-1 




/min 
1 /mill 


N-n 

1 1 




N-n-1 


N-n "' N-n-1 


^n can have the values 0, 1, . 


.. N- 


2 





TABLE in. Cross sections, a, for a particle in direction —1 colliding with particles in the other 
remaining directions of the tetrahedral gas 



Collision 


Cross section 




a 


(-1,-2)^(3,4) 


1 


(-1,-3)^(2,4) 


1 


(-1,-4)^(2,3) 


1 


Othor collisions 
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